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Abstract 

We have searched for direct CP violation in quasi-two-body charmless B decays observed in a 
sample of about 45 million B mesons collected with the BABAR detector at the PEP-II collider. 
We measure the following charge asymmetries in decay: Acp{B^ — * rj'K ) = —0.11 ± 0.11 ± 0.02, 
A C p{B ± -> oj^) = -0.0llg;§? ± 0.03, A C p{B± -> fR*) = -0.05 ± 0.20 ± 0.03, 
A C p{B ± -> <f>K*±) = — 0.43±g;|g ±0.06, and A C p(B°/B° -> <PK*°/K*°) = 0.00 ± 0.27 ± 0.03. 
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1 Introduction 



^CP = |„ |2 - ,„ |2 | n nl — AJL , (2) 



The phenomenon of CP violation has played an important role in understanding fundamental 
physics since its initial discovery in the K meson system in 1964 Q. Recently, a significant CP 
violating asymmetry has been observed in the B meson system [Q. Both effects may be accounted 
for by a non-zero phase in the mixing of two neutral mesons (K° — K° or B° — B°). There is a 
different type of CP violation due to interference among decay amplitudes which differ in both 
weak and strong phases. This "direct" CP violation has been observed recently in kaon decays Q. 
While CP violation effects are small in the kaon system they are anticipated to be larger in B 
decays jjj]. Direct CP violation would be evident in an asymmetry of B decay rates: 

_ T(B^f)-T(B^f) 
CP T(B^f)+T(B^f)- {) 

Rare B meson decays are particularly interesting in searches for direct CP violation because 
they have significant penguin amplitudes. In the Standard Model substantial CP violation in B 
decays could arise from interference of penguin (P) and tree (T) amplitudes Q : 

2 \P\ \T\ sinA</> sinA<5 
\P\ 2 + |T| 2 + 2 \P\ jrj cosA</> cos AS 

where Acj) and AS are the differences in weak and strong phases. The weak phase difference, Acj), 
between the b — ► u tree and b — > s (or b — > d) penguin amplitudes is 7 (or 7 + /?), as in the case 
of the decays B^ — > rj'K (or B — > low ). Thus, Acp is sensitive to the CKM angles 7 and a = 
7r — (7 + 0), where 7 = arg(V^* b ) and = arg(I4d) m the usual phase convention []|, H. However, 
there is large uncertainty in the strong phases, which weakens any quantitative relationship to the 
weak phase angles. 

Even more interesting is the scenario of direct CP violation in the pure penguin modes, such 
as B — > (f)K^*\ In the Standard Model, the expected Aqp is negligible. However, new particles 
in loops, such as charged Higgs or SUSY particles, would provide additional amplitudes with 
different phases. Depending on the model parameters, Acp may be as large as 30% with new 
physics []?]]. Complementary searches for new physics would involve measurements of the time- 
dependent asymmetries in B decays to CP eigenstates, such as (pK^^ and rj 1 K^ L y Comparison 
of the value of sin 20 obtained from these modes with that from charmonium modes can probe 
for new physics participating in penguin loops. In these measurements, direct CP violation in the 
decay becomes highly relevant. 

A search for direct CP violation in B meson decays to ttK, r]'K, and ujtt was performed previ- 
ously by the CLEO experiment || . In this paper we improve the precision of the measurements and 
extend the search for direct CP violation to new modes with data from the BABAR experiment. We 



present measurements of the charge asymmetry in the quasi- two-body charmless B decays ||, 10 1: 
B± j/jfi, B ± w^, B ± -> ^>K ± , B ± -» and B°/B° 4>K*°/K*°. We choose modes 

where the B flavor is tagged by its charge, except for the (f)K*°/K*° final state where the flavor is 
tagged by the charge of the kaon from the K*° — > K + ir~ decay. A measurement from BABAR of 
the ttK charge asymmetry may be found elsewhere [pLlfl . 

2 Detector and Data 



The data were collected with the BABAR detector IS at the PEP-II asymmetric e + e collider [13| 



located at the Stanford Linear Accelerator Center. The results presented in this paper are based on 
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data taken in the 1999-2000 run. An integrated luminosity of 20.7 fb 1 was recorded corresponding 
to 22.7 million BB pairs at the T(4S) resonance ( "on-resonance" ) and 2.6 fb" 1 about 40 MeV below 
this energy ("off-resonance"). 

The asymmetric beam configuration in the laboratory frame provides a boost to the T(4S) 
increasing the momentum range of the S-meson decay products up to 4.3 GeV/c. Charged particles 
are detected and their momenta are measured by a combination of a silicon vertex tracker (SVT) 
consisting of five double-sided layers and a 40-layer central drift chamber (DCH), both operating 
in a 1.5 T solenoidal magnetic field. With the SVT, a position resolution of about 40 pm is 
achieved for the highest momentum charged particles near the interaction point, allowing the precise 
determination of decay vertices. The tracking system covers 92% of the solid angle in the center- 
of-mass system (CM). The track finding efficiency is, on average, (98±1)% for momenta above 
0.2 GeV/c and polar angle greater than 0.5 rad. Photons are detected by a Csl electromagnetic 
calorimeter (EMC), which provides excellent angular and energy resolution with high efficiency for 
energies above 20 MeV Q. 

Charged particle identification is provided by the average energy loss (dE/dx) in the tracking 
devices and by a unique, internally reflecting ring imaging Cherenkov detector (DIRC) covering the 
central region. A Cherenkov angle K-ir separation of better than 4a is achieved for tracks below 
3 GeV/c momentum, decreasing to 2.5<r at the highest momenta in our final states. Electrons are 
identified with the use of the tracking system and the EMC. 

3 Event Selection 

All the selection requirements are identical to those used in the branching fraction measure- 
ments H [10[| . Hadronic events are selected based on track multiplicity and event topology. We fully 
reconstruct B meson candidates from their charged and neutral decay products, where we recover 
the intermediate states ?/ — ► rj^-K^ (r}L n ^) or p°7 (?7p 7 ), ui — ► 7r + 7r~7r , 4> —>■ K + K~ , K* + — > K°ir + 
(K£o) or K + tt° (Kp- + ), K*° K+ir~, p° -> tt+tT , tt° -> 77, tj -» 77, and K° -» K s -> vr+vr-. 
Candidate charged tracks are required to originate from the interaction point, and to have at least 
12 DCH hits and a minimum transverse momentum of 0.1 GeV/c. Looser criteria are applied to 
tracks forming K s candidates to allow for displaced decay vertices. Kaon tracks are distinguished 
from pion and proton tracks via a likelihood ratio that includes, for momenta below 0.7 GeV/c, 
dE/dx information from the SVT and DCH, and, for higher momenta, the Cherenkov angle and 
number of photons as measured by the DIRC. 

We combine pairs of tracks with opposite charge from a common vertex to form K®, <ft, K*°, 
and p° candidates. We further combine a pair of charged tracks with a it or rj candidate to 
select u) or 77' w candidates. The selection of Kg candidates is based on the invariant two-pion 
mass (\M n7T — m#o| < 12 MeV/c 2 ), the angle a between the reconstructed flight and momentum 
directions (cos a > 0.995), and the measured lifetime significance (r/oy > 3). 

We reconstruct 7r° (77) mesons as pairs of photons with a minimum energy deposition of 30 MeV 
(100 MeV). The typical width of the reconstructed 7r° mass is 7 MeV/c 2 . A ±15 MeV/c 2 interval 
is applied to select 7r° candidates. We combine a p° candidate with a photon of energy above 
200 MeV to obtain an rj pi candidate. 

We select <f>, u>, r/, and 77 candidates with the following requirements on the invariant masses 
of their final states (in MeV/c 2 ): 990 < m(K + K~) < 1050, 735 < m(7r + 7r-vr ) < 830, 930 < 
m(777r+7T~) < 990, 900 < m(pj) < 1000, and 490 < 771(77) < 600. The natural widths of the K* 
and p dominate the resolution in the invariant mass spectrum. We require the invariant p mass to 
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be between 500 MeV/c 2 and 995 MeV/c 2 . For K* candidates the Kit invariant mass interval is 
either ±100 or 150 MeV/c 2 §. 

The helicity angle Oh of a <j>, K*, or uj is defined as the angle between the direction of one of 
the two daughters, or the normal to the uj decay plane, and the parent B direction in the resonance 
rest frame. To suppress combinatorial background we restrict the K* + — ► K + ir° helicity angle 
(cos#ff > —0.5). This effectively requires the 7r° momentum to be above 0.35 GeV/c. 

We identify B meson candidates kinematically using two nearly independent variables [12], the 



energy-substituted mass toes = y ( s /2 + Pi • ps) 2 /Ef — p B and AE = {E^Eb — Pi ■ Pb — s/2)/\/s, 
where y/s is the total e + e~ CM energy. The initial-state four-momentum (E^pj) derived from 
the beam kinematics and the four-momentum (Eb,Pb) of the reconstructed B candidate are all 
defined in the laboratory. An alternative to toes is the energy constrained mass toec 5 which is 
obtained from the kinematic fit of the measured candidate four momentum in the T(4S) frame with 
the constraint of the B energy to the beam energy. Both toes an d toec provide almost identical 
background separation, while toec is l ess correlated to AE than is toes- For signal events AE 
peaks at zero and toes and toec at the B mass. 

Monte Carlo (MC) simulation demonstrates that contamination from other B decays is neg- 
ligible. However, charmless hadronic modes suffer from large backgrounds due to random combina- 
tions of tracks produced in the quark-antiquark (qq) continuum. This background is distinguished 
by its jet structure as compared to the spherical decay of the T. To reject continuum background 
we make use of the angle 9t between the thrust axes of the B candidate and the rest of the tracks 
and neutral clusters in the event, calculated in the center-of-mass frame. The distribution of cos 9? 
is sharply peaked near ±1 for combinations drawn from jet-like qq pairs, and nearly uniform for 
the isotropic B meson decays. Thus we require | cos6>t| < 0.9 (0.8 for <pK* + ). We also construct 
a Fisher discriminant |l4|] which combines eleven variables: the angles of the B momentum vector 
and the B two-body decay axis with respect to the beam axis in the T(AS) frame, and a nine bin 
representation of the energy flow about the B decay axis. 



4 Maximum Likelihood Fit 

We use an extended unbinned maximum likelihood (ML) fit to extract signal yields and charge 
asymmetries simultaneously. The extended likelihood for a sample of N events is 

/ M 2 \ N / M 2 \ 

£ = expl-J2J2n ik \ Y[ 1^2 Uik ^ikixf.a) J , (3) 

\ i=l k=l J j=l \i=l k=l ) 

where Vi^(xj\a) describes the probability for candidate event j to belong to category i and flavor 
state k, based on its measured variables Xj, and fixed parameters a that describe the expected 
distributions of these variables in each of the M categories. This probability is non-zero only for 
the right final state flavor (k = 1 for B — » / and k = 2 for B — > /). In the simplest case, 
the probabilities are summed over two categories (M = 2), signal and background. The decays 
with the charged primary daughter fv^ (K or ir^ 1 ) are fit simultaneously with two signal and two 
corresponding background categories (M = 4). These are: B^ — > ry'/i 1 * 1 , uih^, and (fih^. We rewrite 
the event yields riik in each category in terms of the asymmetry A% and the total event yield nf. 
mi = m x (1 + Ai)/2 and = n« x (1 — Ai)/2. The event yields n, and asymmetries Ai in 
each category are obtained by maximizing C (]Tq ]. Statistical errors correspond to unit changes in 
the quantity x 2 = —2ln£ around its minimum value. The significance of non-zero asymmetry is 
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defined by the square root of the change in \ 2 when constraining the asymmetry to zero in the 
likelihood fit. The 90% C.L. limits correspond to a change in \ 2 °f 2.69. When more than one 
channel is measured for the same primary B decay, the channels are combined with \ 2 distributions. 

The probability Vik{xj\ a) for a given event j is the product of independent probability density 
functions (PDFs) in each of the fit input variables xj. These are AE, tties ° r m EC, invariant 
masses of intermediate resonances (r/, uj, <p, K* , and rf), Fisher discriminant, and the <fi and u 
helicity angles for pseudoscalar-vector decays. For the simultaneous fit to the decays with the 
charged primary daughter we include normalized residuals derived from the difference between 
measured and expected DIRC Cherenkov angles for the h^. Additional separation between the two 
final states is provided by AE, where the separation depends on the momentum of the charged 
primary daughter in the laboratory and it is about 45 MeV on average. 

The fixed parameters a describing the PDFs are extracted from signal and background distri- 
butions from MC simulation, on-resonance AE-m^ sidebands, and off-resonance data. The MC 
resolutions are adjusted by comparisons of data and simulation in abundant calibration channels 
with similar kinematics and topology, such as B — > Dtt, Dp with D — > Kir, Kirir. The simulation 
reproduces the event-shape variable distributions found in data. The Cherenkov angle residual 
parameterizations are determined from samples of D° — > K~n + originating from D* decays. 

For the parameterization of the PDFs for AE, mss or m^o, and resonance masses we employ 
Gaussian and Breit-Wigner functions to describe the signal distributions. For the background we 
use low-degree polynomials or, in the case of tties or toec, an empirical phase-space function [17|. 
The background parameterizations for resonance masses also include a resonant component to 
account for resonance production in the continuum. In the B decays into vector-vector states, the 
cos Oh distribution is the result of an a priori unknown superposition of transverse and longitudinal 
polarizations, and thus it is not used for background suppression in the fit. For pseudoscalar- 
vector B decay modes, angular momentum conservation results in a cos 2 Oh distribution for signal. 
The background shape is again separated into contributions from combinatorics and from real 
mesons, both fit by nearly constant low-degree polynomials. The Cherenkov angle residual PDFs 
are Gaussian for both the pion and kaon distributions. The Fisher discriminant is described by an 
asymmetric Gaussian for both signal and background. 



5 Results 

The results of our ML fit analyses are summarized in Table [j]. The signal yields along with branching 
fraction results have been reported earlier H |l0| . In all cases we find significant signal event yields, 
and hence proceed with asymmetry measurements. The dependence of the x 2 on -Acp for each 
decay mode and sub-channel is shown in Fig. |l| and asymmetry measurements are summarized in 
Fig. § We see no significant asymmetries and set 90% C.L. intervals. 

Most of the systematic error contributions relevant to branching fraction analyses cancel for the 
ratio in Eq. |]. Some level of charge asymmetry bias is inevitable as neither the BABAR detector 
nor PEP-II is perfectly charge symmetric. However these effects are mostly very small for the final 
states considered here. Charge biases in the detector and track reconstruction have been studied in 
a sample of more than a billion charged tracks in multi-hadron events. After proton and electron 
rejection we find an asymmetry consistent with zero with an uncertainty of less than 1% for a wide 
range of momenta. Taking into account particle identification requirements, this consistency is still 
better than 2%. The D*^ control sample of kaon and pion tracks is used to estimate systematic 
uncertainties in the asymmetries arising from possible charge biases in the Cherenkov angle residual, 
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Figure 1: Distribution of \ 2 f° r the charge asymmetries Acp in the physically allowed range. Top 
plots: two secondary channels (dashed and dotted lines) are combined to produce a final result 
(solid line); left plot: rj'K^ with rj'^^K^ (dashed) and rj'K^ (dotted); right plot: (pK*^ with 
4>K*± (dashed) and <t>K*± + (dotted). Bottom plots, from left to right: uw ± , <\>K±, and 4>K*°/K*°. 
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Table 1: Results of the ML fits: number of signal events (n s i g ), their charge asymmetry (Acp), 
asymmetry 90% C.L. limits and significance (Sj\). All results include systematic errors. 



Mode 


n s jg 




yu/o U.Li. 




Tj J\ 


49.5i^ ± 1.5 
87.6±^ ± 3.7 


nil i n 11 j_n no 
— U.ll ± U.ll ± U.UZ 

-0.17 ±0.15 ±0.01 

-0.05 ±0.15 ±0.03 


[— U.zo;+U.lwJ 


1 n 
l.U 

1.1 

0.3 


, >-p± 

U>TT 


^ I -D_7 y ± i.y 


n m +0.29 i n r\Q 
— U.U1_q 31 ± U.Uo 


[— U.oU;±U.4bJ 


n n 
U.U 




31 4 + r'n ±23 
ox - -5.9 


—0 05 ± 20 ± 03 


[_n 37- +n 281 


0.2 


<f>K*± + 


4.4t|£±0.4 
7.1±|;| ± 1.2 


— 0.43±g;|g ± 0.06 
— 0.55±g;U ± 0.05 

n O1+0.54 +0.10 
U - Oi -0.43 -0.06 


[-0.88;+0.18] 


1.2 
1.1 
0.6 




20.8ijj;5 ± 1.3 


0.00 ±0.27 ±0.03 


[-0.43;±0.43] 


0.0 



which are found to be less than 1%. 

From these studies we assign a systematic uncertainty of 1% on Acp for all the modes with a 
charged primary daughter: — > rfh , wh , and cph^. For the modes with a K* we account for 
the broader momentum spectrum of the charged daughters and particle identification applied to 
the kaon candidates with a 2% systematic error. All measured background asymmetries and signal 
asymmetries in MC are consistent with zero within statistical uncertainties. 

A different type of uncertainty originates in the ML fit from assumptions about the signal 
and background distributions. We vary the PDF parameters within their respective uncertainties, 
and derive the associated systematic errors in the event yield and its asymmetry. Corresponding 
systematic errors on asymmetry are found to be 2% for r/K and (j)K*°, 3% for unr^ and (pK^, 
and 6% for (pK*^, the latter being dominated by the mode with a 7r°. These systematic errors are 
conservatively estimated and may be improved with higher statistics. 

We combine the correlated (due to selection requirements) and uncorrelated (due to PDF vari- 
ations) systematic errors, and convolute systematic errors into \ 2 distributions in order to obtain 
results with systematics. We also treat the correlated and uncorrelated systematic errors separately 
when we combine the sub-channels. The uncertainties in the final results presented in Table [l] are 
dominated by statistical errors. 
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Figure 2: Results of the direct CP violation search in the B decays into final states rj'K^, uir^, 
<j)K , 4>K*^, and <pK*° / K*° . Points with error bars represent experimental measurements of Acp- 
Solid thick lines delimit the 90% C.L. intervals. For the modes t/K and ujtt^ smaller points with 
error bars show results of the CLEO experiment B . 
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6 Conclusions 



We have searched for direct CP violation in quasi-two-body charmless B decays observed in BABAR 
data. The measured charge asymmetries of the B decays into final states r/K , lott^, <\>K^ , ifrK*^, 
and <pK*° /K*° are summarized in Table |l| and Fig. ^. The 90% C.L. limits rule out a significant 
part of the physical Acp region. 
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